Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 April 2009 (MN 1^1^ style file v2.2) 



Star clusters around recoiled black holes in the Milky Way 
halo 



Ryan M. O'Leary* and Abraham Loebf 

Harvard University, Department of Astronomy, 60 Garden St., Cambridge, MA 02138, USA 



Q 1 April 2009 

o 

(N 

<: 



o 



(N 
> 
(N 

\o 

(N 

O 
oo 
O 



ABSTRACT 

Gravitational wave emission by coalescing black holes (BHs) kicks the remnant 
BH with a typical velocity of hundreds of kms~^. This velocity is sufficiently large to 
remove the remnant BH from a low-mass galaxy but is below the escape velocity from 
the Milky Way (MW) galaxy. If central BHs were common in the galactic building 
blocks that merged to make the MW, then numerous BHs that were kicked out of 
low-mass galaxies should be freely floating in the MW halo today. We use a large 
statistical sample of possible merger tree histories for the MW to estimate the expected 
number of recoiled BH remnants present in the MW halo today. We find that hundreds 
of BHs should remain bound to the MW halo after leaving their parent low-mass 
galaxies. Each BH carries a compact cluster of old stars that populated the core of 
its original host galaxy. Using the time-dependent Fokker-Planck equation, we find 
that present-day clusters are < 1 pc in size, and their central bright regions should be 
unresolved in most existing sky surveys. These compact systems are distinguishable 
from globular clusters by their internal (Keplerian) velocity dispersion greater than 
one hundred km s^^ and their high mass-to-light ratio owing to the central BH. 
An observational discovery of this relic population of star clusters in the MW halo, 
would constrain the formation history of the MW and the dynamics of BH mergers in 
the early Universe. A similar population should exist around other galaxies and may 
potentially be detectable in M31 and M33. 

Key words: galaxies:kinematics and dynamics-galaxies:nuclei-black hole physics- 
gravitational waves-star clusters 



C3 1 INTRODUCTION 

During the final coalescence of two black holes (BHs), 
gravitational waves (GWs) are emitted unisotropically and 
carry away linear momentum, thus kicking the merger BH 
remnant in the oppo site direction (|Peres JiljjJ; iBekensteinI 

1 19731 : iFitchettllTosS '). The resulting kick velocity of typ- 
ically hundreds of kms~^ depends on the mass ratio of 
the BHs as well as t he spin and orientation of the bi- 
nary b efo re coalescence teakcr e t al ]|2006l : ICampanell i et al.l 

l2007bl P: iTichv fc Marronetti ,20071 '). Such kicks can alter 
both the the population of nuclear BHs in galactic bulges 
(iMadau & Quataert^ '2004'; 'Libeskind ct al.' '200^; IVolonteril 

[2007; Schnittman 2007 ; Blccha & Loeb 200 8.) as well as 
the c ore of stars in the bulge itself (jGualandris fc Merritj 

|2008| ). The discovery of a remnant population of recoiled 
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BHs in the present-day Universe can provide a new win- 
dow i nto the merger statistics and spin dis tribution of the 
BHs (|Libeskind et al.ll2006l : IVolonterill2007D . as well as test 
general relativity in the strong regime. 

Recoiling BHs may be detected through their 
unique flaring (Lippai ct al' ' 20081 : IShields fc Bonning||2008l : 
fS chnittman fc K rolik 2008. ) or observed as spatially and 
kinematically offset quasars l|Madau fc QuataertI |2004| ) . If 
the BH binary was surrounded by an accretion disk, the 
ejected remn ant BH woul d carry the disk with it and shine 
as a quasar l|Loebl '2007) . A search for spectral shifts in 
the broad lines of quasars relative to the narrow emission 
line s of their parent gal axies resulted mainly in upper lim- 
its ((Bonning et al."2007') and possibly one suggested candi- 
date (Komossa ct al. 2008) for an ejected quasar. Unfortu- 
nately, these quasars lose their source of gas for accretion , 
and have short Ufetimes (|Loebl l2007l : iBlecha fc Loebl l2008l : 
IVolonteri fc Madau|[200i ). 

Past studies of the observational consequences of GW 
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recoil have generally focused on massive BHs with a mass 
M, > l(f Mq. Such BHs reside in the most massive galaxies, 
and require the greatest and rare st kick velocities in order 
to es cape their parent galaxies (|Schnittman fc Buonannol 
I2OO7I ). But since the kick velocity from BH mergers is in- 
dependent of the total mass of the system, the dynamical 
consequences are most prominent for less massive BHs which 
tend to reside in low-mass galaxie s with shallower potential 
wells and lower esc ape velo cities (iMadau fc Quataertll2004l: 
Volonteri fc Pernal 2005; B aker et al.ll2006l: iLibeskind et al l 



2006l : IVolonterill2007i : iTanaka fc HaimanlbOOSl ). In fact, most 
BHs with M. < lO'^ M© wiU likely be kicked out of their 
paren t galaxies as a result of major mergers (|Baker et al.l 
I2OO6I ). During the hierarchical build-up of the Milky- Way 
(MW) galax;y, mergers of low-mass galaxies were common 
and s hould have r esulted in a population of freely floating 
BHs (|Volonteri fc Pcrna 2005). Since the gravitational po- 
tential of an overdense region does not evol ve dramatic ally 
during the growth of cosmological structure l|Loebll2o"o^ the 
region that eventually collapsed to make the MW was able 
to trap those BHs with the most common kick velocities 
(< 500 km s~^) even before the MW had formed. Each of 
the ejected BHs carries with it a sta r cluster that used to 
popu l ate the core of its pare nt galaxy l|Gualandris fc MerrittI 
I2OO8I : iKomossa et"aLll2008l ). In this paper we examine the 
observational signatures of the ejected star clusters which 
are expected to be floating in the MW halo. The discovery 
of a relic population of the star clusters attached to recoiled 
BHs would provide a new window for cosmology, allowing 
one to constrain the merger history of the MW as well as 
the formation history of the first population of massive BHs. 



Previous studies have looked at populating the 
MW halo with rec oi ling B H s as analysed aga in here 
llMadau fc QuataertI |2004| : IVolonteri fc Pernal l2005l : 
ILibeskind et al.l |2006|). through three body encounters 
( Volonteri fc Perna 2005), as weU as a variety of other pro- 
cesses. Other sources of w andering BHs are t he remnants 
of P opulation HI stars l|lslam et al.l l2003l : IZhao fc Sil^ 
I2OO5I ) as well as the the direct collapse of b aryons into a 
BH (|Bertone et al.l l2005l : iMapelh et al] l2006t ). In both of 
these cases, however, the BH either is presently located 
in the center of a galaxy or is "naked" without any 
stellar companions. Such BHs are likely to only be located 
by their i nteraction with the matter that surroun ds them 
Jlslam et a l. 2004c, b ; Volonteri fc Pcrna 2005; Mil fc Totanil 
I2OO5I : I Bertone et al.l I2OO5I ). Additionally, dwarf galaxies 
may be tidally stripped of their outer stars, leaving behind 
a BH surrounded by a massive cluster of s t ars similar to 
some globular clusters (|Zinnecker et al.lll988l : [Bassino et al.l 
1 19941 ). In contrast, the new source of BHs presented in this 
paper uniquely have a population of bound stars that are 
only a fraction of the mass of the central BH. 



In § [2] we consider the formation history of the M W to 
estimate the number and distribution of ejected BHs that 
may be found in its halo today. We then follow the long term 
evolution of the star clusters that are attached to them in 
§ 131 In § SI we discuss the observational signatures of these 
systems. Finally, § [5] summarizes our main results and their 
implications. 



2 THE MERGER HISTORY OF THE MW AND 
EJECTED BHS IN THE HALO 

In the stan dard cosmological co ntext of hierarchical galaxy 
formation (|Springel et al.l |2008| ) the MW formed from an 
overdense region of the uni verse with a present-day virial 
mass Mvir « 1-2 x 10^^ Mq (|Klvpin et al.ll2002l : iLi fc Whitd 
|2008| ). The threshold for cooling by atomic transitions limits 
the minimum mass o f star forming galaxies at high redshifts 
to Afgai ~ 10* M0 (|Wvithe fc Loebl |2006| ). If the earliest 
population of star-forming galaxies had BHs at their cen- 
ter, then their shallow potential would be insufficient to re- 
tain the recoiling BHs after mergers with a BH mass ratio 
q = (M,,i/M.,2) 0.1. Nevertheless, since most kicks are 
modest (< 500 kms^^), the ejected BHs remain confined to 
the MW. 

We estimate the number of BHs in the MW halo, A^'bh, 
by analyzing a large statistical ensemble of possible merger 
tree histories for the MW using analytic models for the dis- 
tribution of kick velocities and galaxy profiles. In particular, 
we generate > 10'^ Monte-Carlo realizations of the merger 
history of the MW g alaxy with the merger tree code made 
publicly available bv lParkinson et"al] (l2008t). This code fol- 
lows a modified form of the extended Press fc Schechteil 
([19741) formalism and is normalized to match the Millennium 
run n umerical simulation (ISpringel et al ] I2OO5I : ICole et al.l 
l2008h . We assume that the MW has a mass of 1.5 x 10^^ M©, 
and restrict our simulations to formation scenarios in which 
the MW un derg oes no major me rgers {q > 1/3) since a red- 
shift z = l (jPiemand et ai.ll2007^ . 

We determine the prope rties of the halo from its circular 
velocity at the virial radius (|Barkana fc Loebll200ll ). 
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where = IStt^ + 82d - 39d^ d = fi^ - 1, = 
(r2m(l + z)^/(Slm(l -I- z)'^ ^ f^A), evaluated at the merger 
redshift z and we adopt the values for the cosmological den- 
sity parameters, Q.m and ^a, used in the Millennium run 
(|Springel et al.ll2005^ . In each galaxy, we as sume that the 
dark matter follows an NFW radial profile (iNavarro et al.l 
Il996h with a concentration parameter c = 4 out to the virial 
radius, as expected for a newly formed dark-matter halo 
(|Wechsler et al.ll2002l ). In this case, the escape velocity from 
the galaxy's center is iiesc ~ 2.%Vc- We scale the stellar veloc- 
ity dispersion of the bulge in the low-mass galaxies similarly 
to the MW and adopt cr* = fc/2. The BH mass in the cen- 
ter of ea ch galaxy is then dict ated by the observed M,-a* 
relation (|Tremaine et al.ll2002l ) 



km s 



M. = 8.1 X 10^ Mq 
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We assume that the merger of the BHs is efficient, occur- 
ring before a third BH is introduced into the system, and 
that the remnant BH mass is M.^flnai = -I- M,,2. Al- 
though the M,-cr* relation was i nitially determined from 
the analys is of more massive BHs, iBarth et all (|2005l ) and 
I Greene fc H o (2006) have found that it consistently extends 
to active galactic nuclei with BHs with masses ^ 10^ M0, 
even without classical bulges (|Greene et al.|[200^ ). 

For each merger tree, we randomly assign a kick veloc- 
ity from analytic models of the kick velocity distribution of 
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ISchnittman fc Buonannol l|2007l ) assuming that all the BHs 
have the same spin with two different values: a — 0.1 and 
0.9. A BH is ejected from its host galaxy if the kick ve- 
locity is larger than the escape velocity, Vk > Uosc, and 
remains in the larger potential of the MW if the final ve- 
locity is less than the escape velocity of the MW halo, 
Vc} — \^v'^ — Wcsc ^ 500kms~^. Although the MW escape 
velocity evolves by tens of percent as the MW halo assem- 
ble^, we neglect this evolution in our analysis. 

In our realizations of the MW merger tree, there were 
an average of « 1500 galaxy mergers for each run. Of these, 
~ 700 were major mergers with q > 0.1. For BHs with spin 
a — 0.9 (a = 0.1) gravitational wave recoil ejected ~ 570 
(« 440) BHs from their parent galax;y, and A''bh ~ 330 
440) remained within the MW halo. In all of our sim- 
ulations, the average kick velocity, even for the most mas- 
sive BHs ejected, is comparable to or less than the present 
day velocity dispersion of the halo < 200kms~^. Thus, we 
expect that most of the BHs trace the distribution of the 
d ark matter in the halo. This is consistent with the results 
of iLibeskind et akl (|2006l ) who used A'^-body simulations to 
follow the kicked BHs in MW halo. The BHs which were not 
ejected from their host galaxies settle back into the core of 
galaxies owing to dynamical friction, however, the vast ma- 
jority of these galaxies have since merged to form the Milky 
Way. 

In Figure [Tl we plot the cumulative flux distribution of 
the remnant BHs, assuming that the BHs follow the present- 
day NFW profile of the MW with c = 12 and rvir = 200 kpc 
iKlvpin et al.ll2002l : lLrfc Whitell2008l '). We determine the lu- 
minosity from Eq. (|6]) in § [S] and report the flux in units of 
the flux from the Sun if it were at a distance of 1 kpc, /©.kpc, 
and also the apparent bolometric magnitude (top axis). As- 
suming that these clusters lose little mass, we expect that 
nearly all remnant BHs with mass M, > 2 x 10^ M0 in the 
MW halo would be visible to the depth of Sloan Digital Sky 
SurveyE (SDSS). 

In all our models, we find the mass distribution of the 
ejected BHs to scale roughly as dNBu/dM, oc AI^^ (with 
the a = 0.1 case showing a slightly steeper distribution on 
the high mass end of BHs). The most massive BH present in 
the halo is G-StH x 10^ Mq for a = 0.9, and 1.7toi x 
for a — 0.1, where the quoted uncertainty represents the 
one standard deviation. The closest BH that has an appar- 
ent magnitude < 21™, the approximate depth of SDSS, has 
an average mass of ~ 3 x 10"^ Mq and is ~ 1 kpc away. 
For the typical distances and masses of the BHs, dynam- 
ical friction is negligible over the age of the universe. For 
example, the dynamical friction timescale of a ~ 10^ M0 
BH is > 10^^ yr at even ~ 1 kpc from the galactic center 
iBinnev fc Tremaind[l987l ). 

Given the assumptions used here and the uncertainties 
in the lower limit for the BH mass, our results are consis- 
tent with the findings of previous analyse s of BHs in the MW 
halo, even though only a few papers (i.e.. lVolonteri fc Pernal 

icitly account for GW re- 
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Figure 1. The cumulative distribution of ejected star clusters 
in the MW Halo. Plotted is the flux distribution associated with 
BHs masses greater than 10'^ Mq (black), 10* Mq (blue), and 
10^ Mq (red), in our models with BH spin a = 0.1 (solid) and 
a = 0.9 (dashed) lines, plotted in units of the flux of the Sun 
at a distance of Ikpc (/Q,kpc)- The top axis is labelled with the 
apparent bolometric magnitude of the clusters. Nearly all BHs 
with M, ^ 2 X 10^ Mq have apparent magnitudes greater than 
21, the rough magnitude limit of SDSS. The mass distribution 
of the ejected BHs has approximately equal number per logJ\/, 
interval, with dA^Bn/dA/, oc . 



coil. Previous studies that focused on BHs that formed in 
the smallest overdensities in the early universe (- 10*^ Mq 
at 2 > 20) found tha t there may be up t o 10^ - 10'* BHs with 
masses > 10 ^ Mq^ dlslam et all l2004a| : IZhao fc Sil^ l2005l: 



^ The gravitational potential does not evolve during the linear 
growth of perturbations at redshifts 2 ^ 1, an d is only m odestly 
enhanced in the final collapse of galaxy halos llLoebll2006f) . 
^ http://www.sdss.org 



iBertone et al]|2 005': MapcUi ct al.l |2006r i. lVolonteri fc Pernal 
l|2005l ) and ILibeskind et al.. (20061 ) used more conservative 
estimates for the threshold of forming seed BHs, and ac- 
counted for the growth of the BH owing to accretion. They 
concluded that there are more likely only ~ 10^ such BHs 
in the halo. Here, we sidestepped the issue of how and 
when the seed BHs formed, and assumed that all galax- 
ies (with mass > 10* Mq) have a central BH that follows 
the M, — (J* relationship, which consistently extends to the 
lowest mass BHs observed (jGrecno & Ho 200i). With these 
cond itions for the BH mass, o ur res ults are most co n sisten t 
with IVolonteri fc Pernal (j2005h and ILibeskind et all lj2006l l. 
however we do not follow the evolution of BHs with mass 
< 10=* Mq. 



3 STRUCTURE AND EVOLUTION OF 
EJECTED STAR CLUSTERS 

The number of stars that remain bound to the ejected BH 
will depend on the stellar distribution immediately before 
merger and the magnitude of the kick from GWs. The stel- 
lar distribution is determined by the relaxation timescale 
at the radius of infiuence of the BH binary, — GM,/al, 
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l|Bahcall fc Wolj|l976l ) 



3(2na 



2n3/2 



327r2G2mYn.lnA' 

where a* is the stellar velocity dispersion after the galaxies 
merge, m* is the average stellar mass, is the number 
density of stars at r^, and In A ~ ln(M,/M*) is the Coulomb 
logarithm. From the M.-ct* relation (Eq. [2| the relaxation 
timescale can be written as 

5/4 



10= 



/ M. 



VlQS M0 



(4) 



if we assume that = IM©, and that the total mass in 
stars interior of is ~ 2M, . Within a relaxation timescale, 
the stars form a density cusp within the radius of influence 
around the central BH, with n*(r) oc r~", where a = 1.75 
for a population of equal mass stars and a « 1.5 — 2.0 
for varyin g m ass distributions of the stars and compact 
remnants flBahcall fc W olf 1976:, Il977l : iFreitag et all l2006l : 
iHopman fc AlexandeiiboOGb l . If the BHs coalesce on a time 
much shorter than the relaxation timescale of the stars ow- 
ing to a gas rich merger, then the stellar density profile is 
likely much shallower with ami du e to the scattering off 
the i nspiraling binary (M erritt fc Szell .2006: Mcrritt et al, 
I2OO7I ). However, new stars might form out of the gas and 
introduce a new cusp of stars. In either case, however, the 
relaxation timescale of the stellar systems is much shorter 
than a Hubble time, and the present-day distribution of stars 
is different than its initial condition. 

Assuming that the stars follow a power-law density pro- 
file and that the total stellar mass within is 2M, , the 
stellar density profile before the BH ejection is 



n*(r) 



M. 3 



(5) 



The kick is imparted to the BH merger remnant on a 
timescale much shorter than the dynamical time of the star 
cluster. In the frame of the BH, all stars receive kicks with 
a reflex velocity —Vk- Stars with total energies > — m*Wfe/2 
will become unbound to the BH. In the Keplerian potential 
of the BH, this roughly corresponds to stars at r > = 
y/ GMm I v\ = {a* /vkYri. The total number of stars that 
remain bound the the BH is then. 



2M. 



Vk 



(6) 



For a — 1.75 the total number of stars that remain bound 
for a minimally ejected BH with Vk ~ 5.6(t* is then ~ 
4 X 1O^(M./1O^M0). The actual number of stars may be 
less due to the ejection of stars by the B H binary inspiral 
ijMerritt fc SzellH2006l : [Merritt et al.ll2007l ). These stars wiU 
mainly remain within a radius r k ~ (rj-i, /vkY'ri <C 1 pc, an d 
have a a = 4 profile for r > rk (|Komossa fc Merrittll200^ . 
Given the short relaxation timescale of the star cluster and 
the lack of a source of new stars, the star cluster will expand. 

We determine the long term evolution of the ejected 
cluster by numerically solving the time dependent, angu- 
lar momentum averaged, Fokk er-Planck equation for stars 
around a central massive object l|BahcaU fc Wolj|l976l . ll977l ) 



dg{x,T) 



(7) 



where x - 



dr dx 
-E/{mi,a^) is the dimensionless energy, r = t/t. 



is the dimensionless time, g(x,T) — [{2TTa^y'^^^n^^]f{E) is 
the dimensionless distribution function of the stars, Q{x) is 
the rate stars fiow to energies larger th an x, and R(x) is the 
tidal disruption rate of stars by the BH l|Komossa fc Merritt! 
l2008h . In these units. 



Q{x) 



and 



dy[max{x,y)] ^''^ [ g{x 



R{x) = 



dy 



dx 



ln[J,{x)/JLcV 



(9) 



where Jc{x) — GM,a^^{2x)~^^^ is the angular momen- 
tum of a circular orbit at a radius r « ri/{2x), and 
Jlc is the maximum angular momentum for the star to 
be destroyed or consumed by the BH. For a star of ra- 
dius Ri, Jlc = GM,a~^x~J^{2{x + xtd))^^'^, where a;td ~ 
assume that far from the BH the 
stars initially follow a Maxwellian distribution with disper- 
sion a*. In dimensionless units, the distribution function of 
the unbound stars is just g{x) — exp{x) for a; > 0. For a 
power-law distribution function with g{x) oc x^, the density 
profile of the stars will also follow a power-law n* oc r"", 
with a = —p — 3/2. 

We have performed two sets of calculations, first for a 
binary that merges because of stellar relaxation, and then 
for a gas rich merger. For the first case, we evolved the sys- 
tem for one relaxation timescale, at which point the system 
has reached steady state, and then we removed all unbound 
stars by setting g{x) — for a; > 0. The initial kick on the 
system breaks the spherical symmetry assumed here. How- 
ever, after a relaxation timescale we expect the system to 
reach spherical symmetry again. To represent the kick on 
the bound stars, we scale the distribution function of stars 
as g{x) g{x)z^-^/{l + z^'^), where z = x/{vk/at,f. This 
yields an asymptotic density profile wit h n oc for r rk, 
as ex pected immediately after the kick IjKomossa fc Merritt! 
|2008f) . To examine the evolution of the star system if the BHs 
coalesce through a gas rich merger fc allegari erar '2008'l , we 
adopt g{x) oc x~^^^ z'"^ /{l+z'^), normalized by the same imm- 
ber density and relaxation timescale used in Eq. ([Sjl. This 
yields a ~ 1 for r < rk and a ~ 4 for r > rk, as expected fo l- 
lowing a merger (jMerritt fc Szellll2006l : iMerritt et aLlbOOTl ). 
This run results in ~ 10% the number of bound stars as 
with a = 1.75. Its importance, however, is uncertain given 
the short relaxation timescales for these systems. 

Shortly after the GW recoil kicks the BH, the star clus- 
ter expands as a result of two-body relaxation and evolves 
until the relaxation timescale of the cluster becomes compa- 
rable to the age of the system. Figure [2] shows the evolution 
of the projected number of stars interior to a radius r for a 
BH with mass 10^ Mq for both a = 1.75 and a — 1 density 
profiles. Since our simulations indicate little mass loss, the 
condition that the relaxation timescale of the system (Eq.[3} 
be comparable to the age of the Universe implies that the 
cluster should expand by a factor of order {th/tr)^^^ ^ 10. 
Even then, the present day cluster size is less than the initial 
radius of infiuence of the BH, < 1 pc, and much smaller 
than the tidal radius of the system, distinguishing them from 
globular clusters with a similar stellar mass (owing to the 
gravitational binding provided by the central BH). 

In our simulations there is very little mass loss due to 
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Figure 2. The total number of projected stars interior to r, Np{< r), for M, = IO^Mq and a = 1.75 (left) and a = 1 (right). The 
solid line corresponds to the cluster immediately after being ejected from its parent galaxy with vj. = 5. Scr^ , normalized to have the same 
number density at the radius of influence. The alternating dashed and dotted lines correspond to the projected number of stars after 
every lOtr ~ 650 Myr. Immediately after ejection, the cluster rapidly expands until its relaxation time becomes comparable to the age 
of the Universe, with very little mass loss, n For a 10^ Mq BH, the circular velocity of the stars is ~ 66 kms~^ (r/0.1 pc)"^'''^. 



the tidal disruption of stars or the ejection of stars from the 
system. Since the cluster expands so rapidly, the rate of tidal 
disruption of stars rapidly decreases, and is much lower than 
the rate associated with a comparable BH in the nucleus of 
a galaxy. 

Our simulations only treat the evolution of the cluster 
due to many small-angle scatterings (two-body relaxation), 
and does not include the strong scattering that might oc- 
cur in a dense stellar environment. Such scatterings would 
launch stars on eccentric orbits that take them outside of 
the cluster to about the tidal radius, thus enhancing the 
mass loss from the system. The existence of a halo of such 
stars provides an additional signature unique to the com- 
pact star clusters associated with recoiled BHs. We also ig- 
nore the effects of resonant relaxaation, which can deplete 
the num ber of stars both before and after the ejection of o f 
the BH jRauch fc Tremaind [l99^ : iRauch fc Ingallj Il998h . 
Although resonant relaxation has been approximately ac- 
counted for in Eg .[7] by avera ging over angular momentum 
l|Hopman fc Alexandei|[2006al ) , a full multidimensional anal- 
ysis accounting for angular momemtum is best suited for 
these clusters. 



4 OBSERVING EJECTED SYSTEMS 

Size and Structure: Based on the results in § [2] we esti- 
mate that there should be ^ 300 recoiled BHs in the MW 
halo. Equation (|6]) implies that these BHs are surrounded 
by ~ 4 X 10^(M./10^m*) stars with a half mass radius, 
< 1 pc. At the typical distance to the most massive re- 
coiled BH, ~ 100 kpc, the cluster will have a typical angular 
size < 1 arcsec, below the resolution limit of most existing 



sky surveys. Thus, the cluster will appear as a point source 
and might be confused as a single foreground star with un- 
usual colors. The closest clusters are ~ 1 kpc away and may 
be a few arcmin in size. However, these clusters surround 
the smallest BHs (~ 10'' Mq) with the fewest number of 
stars, and their present day distribution is uncertain owing 
to strong encounters between the stars in the cluster as well 
as stars in the disk. To date, no star cluster has been iden- 
tified in the MW with such a small size; the smallest sta r 
cluster discovered has a size of ~ 3 pc (|Koposov et al.ll2007l ). 
Clusters at distances of tens of kpc could in principle be dis- 
tinguished from point sources by comparing their extended 
image to the point-spread-function of the telescope. In com- 
parison, the stellar mass of globular clusters and the tidally 
stripped cores of dwarf galaxies completely outweights the 
total mass of any central BH in such systems. By measur- 
ing the total mass of stars one can approximately determine 
how large the kick velocity may have been. 
Color: A recoiled star cluster can be distinguished from sin- 
gle foreground stars by its anomalous color. The cluster is 
likely old, ~ 1 — 10 Gyr, and to first order should have colors 
similar to globular clusters of a similar age and metallicity. 
Speetra: The velocity dispersion profile of the ejected star 
clusters should be Keplerian. The unusually large width of 
their spectral lines (greather than one hundred km s~^) can 
be used to distinguish these clusters from foreground stars, 
which possess much lower rotational velocities. 

The compact star clusters around ejected BHs are dis- 
tinguishable from cores of globular clusters by their internal 
(Keplerian) velocity dispersion of order one hundred km s~^ 
and their high mass-to-light ratio owing to the central BH. 
Proper Motion: The most massive systems that reside near 
the virial radius of the galaxy and have a relatively small 
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proper motion. Nevertheless, future surveys like GAIAH, 
Pan-STARRfl and LSSlfl, could aim to measure the 
proper motion of such clusters, and be used in conjunction 
with expected cluster colors to distinguish the clusters from 
foreground stars. The expected proper motion of these clus- 
ters is ~ 10~^ arcsec yr~^. 

The Local Group: Given the constrained distances to M31 
and M33, color-magnitude diagrams could be used to iden- 
tify candidate star clusters around ejected BHs, and spectro- 
scopic follow up could reveal the high dispersions expected 
from the stars very near the BH. At the distance of M31, 
a cluster with rh ^ 0.5 pc would have an angular size of 
~ 0.1 arcsec, and could in principle be resolved by HS10, 
JWSlllI, or ground-based telescopes with adaptive optics. 
Gas Accretion: A faint radio glow from the central black 
hole may also be detectable as it accretes material from 
stella r winds in the surrounding cluster of stars (|Loebl 
bOO^ I. Rarely, these BHs may accrete gas from the inter- 
stellar medium or mol ecular cloud s as they pass through 



the disk of the galaxv llFuiita et al.llT9 98: Isla m et al.ll2003l: 



Mii &: Totaiiil l2005l : IVolonteri fc Pernal 120051 : IMapelli et al.1 
20061 ). however, for the number of BHs ex pected in the MW, 
such a scenario is unlikely to be observed (jVolonteri fc Pernal 
l2005l ). 

Dark Matter Annihilation: The recoiled BHs may be sur- 
rounded by a dense cusp of dark matter, which was adia- 
batically compressed by the baryonic condensation and BH 
growth in its parent low-mass galaxy. These high density 
cusps can dramatically increase the rate of dark matter anni- 
hilation, and may be visib le through the relat i vistic by prod- 
ucts of the annihilation l|Zhao fc "Sii^ l2005l : iBertone et al.l 
l2005l : ISpringel et alllioog) . 



5 SUMMARY AND DISCUSSION 

Based on a large statistical ensemble of merger tree histories 
for the MW, we have found that hundreds of GW recoiled 
BHs should reside within the MW halo today. The BHs have 
a mass M, > 10^ Mq , and should be surrounded by a clus- 
ter of stars that were tightly bound to the BH when it was 
ejected. The most massive BH weighs ~ 1 — 6 x 10^ Mq, 
and is surrounded by a compact star cluster < 1 pc in size 
with a stellar velocity dispersion of hundreds of km s~^. 
High-resolution adaptive optics imaging along with spec- 
troscopy of the star clusters can constrain the distance and 
mass of the BH in the center of the clust er based on orbits 
of individual stars for t he closest clusters (|Ghez et al.ll2005l : 
lEisenhauer et al.ll2005l ). 

The number of recoiled BHs in the MW is most sensitive 
to the fraction of low mass galaxies that harbored BHs in 
their centers, as well as to the merger history of such galax- 
ies during the formation of the MW. We emphasize that 
most of these galactic "building blocks" were made at high 
redshifts and were disrupted during the assembly process of 
the MW; hence, they likely had different properties than the 



http:/ /www.esa.int/esaSC/ 
http: / /pan-Starrs. ifa.hawaii.edu/ 
http:/ /www. lsst.org/ 
http:/ /www. stsci.edu/hst/ 
http: / /www. stsci.edu/jwst/ 



remaining dwarf satellites of the MW today. The future dis- 
covery of the population of ejected star clusters will provide 
a unique probe of the early history of the Milky Way, as 
well as the distribution and evolution of low mass BHs. It 
would also open a new window to exploring the low-mass end 
of the popula t ion of nuclear BHs in high-redshifts galaxies. 
iMerritt et al.l (|2008l ). whose work we learned of after sub- 
mission of this manuscript, arrived at a similar conclusion 
independently, with a slightly different application to the 
nearby Virgo cluster. 
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